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We present an analytical solntion of the single photon quantum feedback in a cavity quantum 
electrodynamics system based on a half cavity set-up coupled to a structured continuum. The 
exact analytical expression we obtain allows us to discuss in detail under which conditions a single 
emitter-cavity system, which is initially in the weak coupling regime, can be driven into the strong 
coupling regime via the proposed quantum feedback mechanism [Carmele et al, Phys.Rev.Lett. 110, 

013601]. Our results reveal that the feedback induced oscillations rely on a well-dehned relationship 
between the delay time and the atom-light coupling strength of the emitter. At these specific values 
the leakage into the continuum is prevented by a destructive interference effect, which pushes the 
emitter to the strong coupling limit. 
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I. INTRODUCTION 


The basic phenomenon at the heart of any quantum 
information processing network is the coherent exchange 
of photonic and atomic excitations by means of a single 
emitter in a microcavity. Advances in the design and fab¬ 
rication of microcavities allow very high quality factors 
and have enabled multiple studies of cavity quantum elec¬ 
trodynamics (cQED) in the strong coupling limit M- 
Several applications are proposed and already realized, 
such as single-photon transistors, two-photon gateways, 
parametric downconversion as well as the generation and 
detection of individual microwave photons 0-[I3- Fur¬ 
thermore, several quantum gate proposals rely on a nat¬ 
ural quantum interface between flying qubits (photons) 
and stationary qubits (e.g. atoms). Here, the photons 
allow for secure quantum communication over long dis¬ 
tances, whereas atoms can be used for the manipulation 
and storage of quantum information 0 - 
Applying cQED techniques require that a single¬ 
atom/single-photon coupling exceeds any photon loss 
and radiative decay processes, such as spontaneous emis¬ 
sion or photon leakage. So, besides technological progress 
to increase the quality factor of the cavities, a promising 
alternative is to identify strategies to control and exploit 
potentially advantageous properties of the environment 
coupling, which go beyond the conventional effects of the 
environment like dissipation and undesired information 
loss. 

A possible mechanism to stabilize qubits and desired 
quantum states is quantum feedback based on the re¬ 
peated action of a sensor-controller-actuator loop. In 
such a case, a quantum system is driven to a target 
state via the external control [ll|, EIIj such that con¬ 
tinuous measurements allow to stabilize the target state, 
e.g. by a modification of the pumping strength. In addi¬ 
tion to these extrinsic control set-ups, experiments start 


to explore a variety of intrinsic, delayed feedback control 
schemes, e. g., by using an external mirror in front of a 
nanocavity |I3| . Intrinsic quantum feedback is not based 
on a continuous measurement process, but controls the 
quantum state by shaping the environment appropriately 

Here, we discuss the photon leakage mechanism shaped 
by an external mirror and show analytically, how the ini¬ 
tial weak atom-cavity coupling is driven into the strong 
coupling regime, following Ref. [2^. Our proposed con¬ 
trol scheme has potential applications for quantum error 
correctio n 1211 , quantum gate purifying |22j| , or quantum 
feedback |1^ . 
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Figure 1: Implementation of an intrinsic quantum feedback 
mechanism via a quasi-continuum, realized by a photonic 
crystal waveguide with length L, which is supposed to be 
considerably larger than the cavity length L'. The waveguide 
is a half-cavity and allows to exchange cavity photons with 
waveguide photons due to the photon leakage G(k). The pho¬ 
tons inside the cavity interact with a single emitter (coupling 
strength 7 ). 
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II. MODEL 


The system consists of a microcavity system of length 
L' with a two-level emitter coupled to a single cavity 
mode (see Fig.[T]). Furthermore, the cavity exhibits pho¬ 
ton loss due to its coupling to external modes. An ex¬ 
ternal mirror, placed in a distance of L, introduces a 
boundary condition to the external mode structure and 
causes a feedback of lost cavity photons into the cav¬ 
ity. We assume that the microcavity length L' L is 
very short in comparison to L to allow a single mode de¬ 
scription for the emitter-cavity interaction. This kind of 
quantum self-feedback can be realized via a shaped mode 
continuum in a photonic waveguide. Due to the finite 
cavity-mirror distance L and the quasi-continuous mode 
structure of the semi-infinite lead, a delay mechanism is 
introduced into the system at r = ^ with cq being the 
speed of light in vacuum. To describe the corresponding 
physics we work with the following Hamiltonian within 
the rotating-wave and dipole approximation : 

H/h = —7 ((T“a^-I-( 1 ) 
- y dfc G(fc, t) + G* (fc, t) dla, 

where a rotating frame is chosen in correspondence to the 
free energy contribution of the Hamiltonian. The emitter 
is described via the Pauli matrices with being the 
raising and lowering operators of the two level system, 
respectively. In the following, the atomic energy is as¬ 
sumed to be in resonance with the single cavity mode. A 
photon annihilation (creation) in the cavity is described 
with the bosonic operator 01(0) and 7 is the coupling 
between the two-level system and the cavity mode. The 
coupling strength between the emitter and the field mode 
is assumed to be of the order of 7 = 50 /reV [ 1 , [H, [2^ . 
The cavity photons interact with the external modes 
in front of the mirror via the tunnel Hamiltonian coupling 
elements G{k,t). Due to the rotating frame and the in¬ 
terference with the back-reflected signal from the mirror, 
these coupling elements depend both on time t and on 
the wavenumber k, resulting in the following expression 
for G{k,t) = Go sin(feL) exp[f(a;o — ujk)t] with Go being 
the bare tunnel coupling strength, loq and uik stand for 
the frequencies of a single cavity mode and half-cavity 
modes, respectively. As we will see below, this specific 
form of G(fc, t) will determine the nature of the feedback 
on the cavity. 


A. Single photon limit 


If no other loss channels or pump mechanism are in¬ 
troduced, the system dynamics described by the Hamil¬ 
tonian Eq. © can be solved in the Schrodinger picture, 
following Ref. [2^ [l^. In the single photon limit, the 
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Figure 2 : Excited state density \ce{t)f of the two-level system 
(black dashed line) and the photon density inside the cavity 
|Cg(t)p (orange solid line) in the quasi-continuum model. The 
quantum feedback mechanism (ft/y = 2) induces a regular 
oscillation pattern at multiples of t = 2n/y. 
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total wave function reads: 

Ivk) = ce|e,0,{0}) + cg|5,l,{0}) (2) 

-k y dfc Cg,fe | 5 , 0 ,{fc}), 

where |e, 0, {0}) denotes the excited state of the two-level 
system with the cavity and the waveguide being in the 
vacuum state, \g, 1, {0}) stands for a single photon resid¬ 
ing in the cavity and the two-level system as well as the 
radiation field in the waveguide being in the ground state. 
Finally, \g, 0 , {fc}) describes the ground state of the two- 
level system with exactly one photon in the waveguide 
of mode k. The variables Ce, Cg, Cg^k denote the corre¬ 
sponding amplitudes of the above three different states. 
Applying the Schrodinger equation, we arrive at the fol¬ 
lowing set of linear partial differential equations: 


^tC-e 

= Cg 

( 3 ) 

dtCg 

= i ^ Ce + i J dk G{k, t) Cg^k 

( 4 ) 


= i G*{k,t) Cg. 

( 5 ) 


First, we numerically solve this coupled set of differen¬ 
tial equations assuming that initially at to = 0 the TLS 
is in the excited state, Ce(to) = 1, and there are nei¬ 
ther photons inside the cavity, Cg(to) = 0, nor in the 
external region, Cg_fc(to) = 0 . To introduce a delay time 
corresponding to r = 2L/co = STr/y, we choose a mirror 
resonator distance L = Trco/y. 

The results for the dynamics of the excited state and 
photon density are shown in Fig. [21 In the time inter¬ 
val [0,t] we find the conventional exponential decay as 
described by the Wigner-Weisskopf model in the weak 
coupling limit. After the first round trip r = 2L/co, the 
photon density and after a small delay, also the excited 
state density, are driven by the quantum feedback. In this 
time interval [r, 2 t] , the amplitude of the photon density 
is smaller than the amplitude of the excited case, since 
the damping mechanism acts only on the photons inside 
the cavity. However, for longer times, the asymmetry be¬ 
tween the amplitudes of the excited state and the photon 
density vanishes, so that the system sets into a state of co¬ 
herent Rabi oscillations characterized by approximately 
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equal maxima of both densities (see the asymptotic dy¬ 
namics for t/r > 8 in Fig. [5]). In this long-time limit, the 
amplitude for the cavity photon population stabilizes at 
around 15% of the maximum photon population in the 
first time interval ([0, r]). This remarkable effect has been 
reported in another publication and will now be an¬ 
alyzed analytically and thereby explained in more detail. 
In particular, we will focus on the following two specific 
questions: (i) How sensitive is this effect on the chosen 
parameters, in particular on the choice of the time de¬ 
lay? (ii) Can the oscillation amplitude be increased or is 
there an intrinsic limit? To answer these questions, we 
will first derive a simplified picture of the dynamics by 
solving Eqs. ([S])-® in the Markovian limit. 



Figure 3: Comparison between numerical calculation and the 
Laplace transformed analytical solution valid until t = 3 t 
with fi :/7 = 2 . 


B. Analytical Quantum Feedback 


mechanism: 


The initial decay and the subsequent oscillations ob¬ 
servable in Fig. [2] indicates that the underlying physi¬ 
cal processes which govern this system consist both of 
a typical (Markovian) cavity loss as well as of a (non- 
Markovian) memory kernel with significant contributions 
around multiples of the delay time r. Assuming that the 
rotating wave approximation and the quasi-continuum 
assumption hold, the coupling to the external modes can 
be eliminated from the problem. To achieve this, Eq. m 
is integrated formally and inserted into Eq. o resulting 
in the following expression: 

dtCg = i ^ Ce — K Cg + K Cg Q{t — t) ( 6 ) 

with K = 7rG'o/(2co). This reduced expression has the 
advantage of being easily solvable numerically and of be¬ 
ing amenable to an analytical solution through a Laplace 
transformation. With the initial conditions, that neither 
cavity nor continuum photons are present in the begin¬ 
ning, i.e. Ce(0) = 1, the equations read after Laplace 
transformation: 

s Ce(s) = 1 -Li 7 Cg{s) (7) 

scg(s) = f 7Ce(s)-KCg(s)-L«:Cg(s) (8) 

where s is the complex frequency parameter of the 
Laplace transformation. As can be seen from Eq. ®, 
the solution consists of a dynamical component without 
the mirror induced feedback t < t and one with the feed¬ 
back for t > r. 

First, we now derive a solution for the photon-assisted 
ground state for t < t, which is the cavity damped 
Jaynes-Cummings model [^ : 

c (A =_ n _=_ ii _ 

^ S 2+72 + KS (s + k/2)2+72-«;2/4- 

This leads directly to the damped Jaynes-Cummings so¬ 
lutions in the time domain as expected for times t < r. 


Cg{t) = i 


Sin 


y/l - (k/ 27)2 7t 


p-k/ 2 i 


v'1-W27)^ 


( 10 ) 


Note, due to the cavity damping not only the amplitude 
is reduced but also the Rabi oscillation frequency is re¬ 
duced by a factor of a/W^(k/ 27 ^. The cavity loss leads 
inevitably to an effectively reduced value for the coupling 
strength, and as a result, the frequency of damped Rabi 
oscillations decreases. As we will see below, this restric¬ 
tion is lifted if a feedback mechanism is present. 

Now, we solve the dynamics for times t > t. This leads 
to an additional term in the denominator. By using a 
geometric series expansion, i.e. (1 — g)“^ = 9 " 

g < 1 and n —>■ 00 , Eq. m can be written as: 


Cg(s) = i 7 


1 K s exp[—(s—'ia;o)T] 

(s+k,/2)2+72 — 


-1 


= % 


(s + k/2)2 + 72-k2/4 


( 11 ) 


[ K s exp[—(s—2 Ci;o)t] 
(s+k/2)^+7^-k=^/4 


(s -f k/2)2 -I- 72 — k2/4 


Due to the linearity of the Laplace transformation, the so¬ 
lution in the time domain can be obtained via the method 
of partial fraction expansion and the convolution prop¬ 
erty. However, the expression is very lengthy and must 
be calculated for every time interval [nr, (n -|- l)r], sep¬ 
arately. Since we are interested in the weak coupling 
regime, we can choose the parameter to be: 7 = k/2 to 
simplify the expression into: 

^ ' ' n 

Using now the binomial series and Laplace transforma¬ 
tion: n!/(s — a)”+^ ^ t" exp[at], we get an expression in 
the time domain: 




(s -I- k/2)2 


( 12 ) 
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Cg{t) 
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-nT)Y, 

fe =0 


(-l)fe [K/ 2 (t-nr)]"+^+^ 
k\{n — ky. (n + 1 + A:)! 


(13) 
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Figure 4: Comparison between numerics (orange curve) and 
long time solution determined by the residuum contribution 
only (black curve) for an initially excited TLS Ce(0) = 1 and 
k /7 = 2. After several r, the analytical long time solution 
and the numerics coincide. 


In Fig. 131 the numerical solution of Eq. (jS]) coupled with 
Eq. ([ 6 |) is compared with the analytical solution Eq. CSD 
for the time interval [0,3t]. The excellent agreement 
found between these two solutions confirms the validity 
of our calculations. For longer times, more terms from 
the series expansion uni) should be included, which is a 
straightforward procedure. As a next step, we derive the 
long time behaviour using the residuum method. 


C. Long time solution 

The long time dynamics of the coupled system is di¬ 
rectly related to the singularities in the contour integral 
of the Laplace transformed function [2^ . To demonstrate 
this explicitly, we need to find the singularities of the 
photon-assisted ground state amplitude: 

I'-y 

9^ ) g2 _|_ ^2 _|_ _ ^gg-(s-iaJo)'r ^ ^ 

The singularities are found by setting the denominator 
to zero. We assume a pure oscillation behavior in the 
long time limit, i.e., where s is purely imaginary. We set 
s = ±17 and get: 

- 7^ -b 7^ ± ijK (1 - = 0, (15) 

from which immediately follows 

gi(<voT7)r = 1. (16) 

We now need to find a delay time r in a way, that for ±* 7 , 
this equation is valid. As it turns out, the corresponding 
singularity condition can be matched for the following 


two cases: 


(f) exp(i 7 T) = exp(iajo'r) = 1 (17) 

(ii) exp{—i^T) = exp(ia;oT) = 1 . (18) 

In order to satisfy these conditions we observe that we 
can freely choose the delay time with respect to the cou¬ 
pling strength by adjusting the length L between cav¬ 
ity and mirror accordingly. For instance, if we choose 
7 T = 27rm and at the same time tune the resonance 
frequency such that ujqt = 27r/, where l,m are integer 
numbers, then the conditions (i) and {ii) are satisfied. 
(Note that there are also three other obvious constraints 
on r and 7 to satisfy the conditions (i) or (ii) but they 
are not discussed below.) As a result we achieve a purely 
coherent asymptotic solution with a minimum of dephas¬ 
ing and a maximum amplitude, corresponding to the fact 
that the pole does not contain any decaying term. In¬ 
deed, we derive the following expression for the asymp¬ 
totic behaviour of the photon density inside the cavity 

c'yKt) = ^ /ds Cg(s) e"* « Res [cg(s)e®*](19) 
(s ± fy) exp[s 2 n 7 r/ 7 ] 17 

^ s^±i'y (s ± i^){s — iy) ± Ks (1 — exp[—s2m7r/7]) 


Using now L’Hopital’s rule and taking the limits s 
± 17 , the solution for Cg{t) reads 


Cg{t) 


i sinjyt] 

1 ± Kmir/j 


( 20 ) 


In Fig. |4j the numerical solution and the analytical long 
time solution is plotted for r = 2tt/j (i.e. when m = 1 ) 
up to several r. The agreement is excellent with the long 
time solution accurately recovering the amplitude and 
the oscillation frequency of the numerical solution. In¬ 
terestingly, for K = 0 , we recover the Jaynes-Cummings 
solution as in Eq. m with K = 0. In contrast to Eq. 
we see, however, that the cavity loss does not modify 
the frequency of vacuum Rabi oscillations, which is now 
equal to the coupling strength 7 , and only damps the 
corresponding amplitude. In this context, we discover a 
maximum amplitude for the feedback effect via this pro¬ 
posed mechanism. It is seen that Eg. (1201) yields the fol¬ 
lowing amplitude of the quantum feedback induced Rabi 
oscillations for r = 271 / 7 : 1/(1 ± KTr/y)^. Therefore with 
K = 27 , the maximum amplitude is approximately 0.02 
in correspondence with Fig.lH which is 15% of the max¬ 
imum photon amplitude in the first time interval [0,t]. 
With these results at hand we can now also answer the 

















questions raised above: 

(I) The effect of stabilized Rabi oscillations in the long 
time limit depends strongly on the chosen time delay r, 
which has to be chosen such to satisfy one of the condi¬ 
tions (El), dHI) that lead to asymptotically undamped 
Rabi oscillations. Furthermore, the factor expliwoT"] 
plays a crucial role to decide whether quantum feedback 
leads to a stabilized Rabi oscillation or to a damped feed¬ 
back situation. However, the effect depends only quan¬ 
titatively (rather than qualitatively) on the cavity loss n 
and coupling strength 7 , besides the obvious restriction 
that both of them are unequal to zero. 

(II) For a given ratio between the coupling strength and 
the cavity loss, x = re/y, there is a maximum amplitude 
which is given for the above case by (1 -I- xtt)”^. 



Figure 5: Fast stabilization of Rabi oscillations of the cavity 
photon number in the strong coupling limit 7 = 20 k after 
only one roundtrip r, when the feedback from the waveguide 
is present (orange curve). Initially, the excitation is stored in 
the TLS (ce(0) = 1). 


D. Photon-Path Representation 


To give an intuitive explanation for this effect of re¬ 
covered Rabi oscillations in the weak coupling limit, we 
visualize the resulting cavity dynamics in the framework 
of the photon-path representation MM- For this pur¬ 
pose we rewrite the system of equations of motion ([7|), 
(IH) in the Laplace domain as: 




with the scattering matrix: 


L = 




( 21 ) 


( 22 ) 


Assuming that there is an inverse to that matrix and 
using the Neumann expansion, we get for ||L|| < 1: 


( Ce(s) 

VC9(S) 



(23) 



Now, the dynamics is written in a very complicated man¬ 
ner but in a way that the photon path (represented by 
scattering processes due to L) becomes visible. With this 
expansion, one can represent the dynamics as a series of 
single scattering events by multiple application of the ma¬ 
trix, which swap the excitation from Ce to Cg and includes 
the cavity loss and the gain from the feedback. This be¬ 
comes especially apparent, when writing down the single 
terms of the Laplace transform and then transforming 
them back into the time domain. In particular for the 
ground state probability such an expansion reads (terms 
up to are kept only): 




(iy) {i'y)K 2 , (*7) 


2 ! 


r + 


3! 


+ 


(»7)'^ 

3! 








(2i7)k^ 


{t - T)'^9{t - t)+... (24) 


From the structure of this expansion follows that un¬ 
damped Rabi oscillations can be viewed as a result of 
an interference between incoming and outgoing photonic 
paths provided that r = 27r/7. In other words, the 
strong coupling feature is produced by a destructive in¬ 
terference effect of the photon paths at the point within 
the waveguide, where the tunneling event between cavity 
and waveguide takes place. This expansion explains fur¬ 
thermore that it takes a minimum time for this effect to 
unfold, since at least two dissipatively interacting waves 
need to be in the waveguide. 


E. Strong coupling limit 


To complete the picture, we investigate the proposed 
feedback mechanism via a quasi-continuum in the strong 
coupling limit. In Fig. [SJ the dynamics with and without 
feedback is plotted for a a coupling strength 7 = 20k. 
Clearly, Rabi oscillations are visible with and without 
feedback. If no feedback is present, however, the am¬ 
plitude of the Rabi oscillations are damped fast without 
changing the frequency. With a feedback and a chosen 
delay time of r = on the other hand, the ampli¬ 

tude loss is stopped at very early times already: Already 
after one roundtrip the amplitude stays constant for all 
times, if no other loss mechanism inside or between mir¬ 
ror and cavity is present. We explain the acceleration of 
the stabilization feature by the fact, that for the strong 
coupling regime the in- and outgoing photons already in¬ 
terfere efficiently after one roundtrip. In contrast, in the 
weak coupling limit the in- and out tunneling does not 
overlap for the first three roundtrips and as a consequence 
inteference takes place at longer times only. If we choose 
a larger roundtrip time r, also in the strong coupling 
regime, a higher number of roundtrips (nr) is necessary 
to reach the point of stabilized Rabi oscillations. 
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III. CONCLUSION AND OUTLOOK 

We have discussed an approach to stabilize single¬ 
emitter cQED via a quantum feedback mechanism in¬ 
duced by an external mirror. Our analytical solution 
shows that depending on the chosen parameters an in¬ 
trinsic limit of the feedback effects exists. For a system 
initially in the weak coupling regime (before the feed¬ 
back modifies the system dynamics) we demonstrate that 
the quantum feedback can at most recover approximately 
15% of the maximum cavity occupancy in the first time 
interval. Our analytic calculations demonstrate further¬ 
more that the quantum feedback induced Rabi oscilla¬ 
tions are indeed coherent and follow a typical differen¬ 
tial delay equation with an appropriate inhomogeneity 


to drive the system into the strong coupling regime. Our 
results extend the set of exact analytical solutions in the 
field of coherent atom cQED and form a starting point 
to establish a framework for a theoretical description of 
coherent quantum feedback. 
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